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Abstract Recent measurements of nucleon resonance transition form factors with CLAS at Jefferson Lab are dis- 
cussed. The new data confirm the assertion of the symmetric constituent quark model of the Roper as the first radial 
excitation of the nucleon. The data on high Q 2 n7r + production better constrain the branching ratios /3ntv and pNr, 
for the Sn (1535). For the first time, the longitudinal transition amplitude to the Sn (1535) was extracted from the 
n7T + data. Also, new results on the transition amplitudes for the Di3(1520) resonance are presented showing a rapid 
transition from helicity 3/2 dominance seen at the real photon point to helicty 1/2 dominance at higher Q 2 . I also 
discuss aspects of the search for new excited nucleon states in strangeness photoproduction. 
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1 Introduction 

Electroexcitation of nucleon resonances has long been 
recognized as a sensitive tool in the exploration of the 
complex nucleon structure at varying distances scales. 
Resonances play an important role in fully understand- 
ing the spin structure of the nucleon. More than 80% 
of the helicity-dependent integrated total photoabsorp- 
tion cross section difference (GDH integral) is the result 
of the iVA(1232) transition^ ^, and at a photon virtu- 
ally Q 2 — 1 GeV 2 more than 50% of the first moment 
rf (Q 2 ) = J Q g 1 (x,Q 2 )dx of the spin structure function g 1 
for the proton is due to contributions of the resonance 
region at W < 2 GeV ^ , and is crucial for describing the 
entire Q 2 range of rf (Q 2 ) and T p - n (Q 2 ) for the proton 
and proton- neutron difference respectively Nu- 
cleon resonances may also be responsible for significant 
portion of the orbital angular momentum contributions 
to the nucleon spin due to the strong excitation of reso- 
nances with quark orbital angular momentum AL Z = 1 
and AL Z — 2 in the 2nd and 3rd resonance region. 

Studies of nucleon resonances are of interest in their 
own rights. Studying the systematics of the excitation 
spectrum is key to understanding the symmetry under- 
lying the nucleon matter. Electroexcitation of resonances 
allows us to probe the internal structure of the excited 



state knowing the structure of the ground state. The 
most comprehensive predictions of the resonance excita- 
tion spectrum come from the various implementation of 
the symmetric constituent quark model based on broken 
SU(6) symmetry ^ . Other models predict a different 
excitation spectrum, e.g. through a diquark-quark pic- 
ture, or through dynamical baryon-meson interactions. 
The different resonance models not only predict differ- 
ent excitation spectra but also different Q 2 dependence 
of transition form factors. Mapping out the transition 
form factors will tell us a great deal about the underly- 
ing quark or hadronic structure. 

The CLAS detector is the first full acceptance instru- 
ment with sufficient resolution to measure exclusive elec- 
troproduction of mesons with the goal of studying the 
excitation of nucleon resonances in detail. With the GeV 
electron beam energy available at JLab, the entire reso- 
nance mass region and a wide range in the photon virtu- 
ality Q 2 can be studied, and many meson final states are 
measured simultaneously^ 1 . In this talk I discuss recent 
results from the electroproduction of single pions to de- 
termine the transition amplitudes to several well-known 
states. I will also discuss selected aspects of the search 
for new excited baryon states. 
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2 The 7VA(1232) transition 

An interesting aspect of nucleon structure at low en- 
ergies is a possible quadrupole deformation of the nu- 
cleon's lowest excited state, the A(1232). Such a de- 
formation would be evident in non-zero values of the 
quadrupole transition amplitude E 1+ from the nucleon to 
the A(1232) EI . In models with 577(6) spherical symme- 
try, the NA transition is simply due to a magnetic dipole 
M 1+ mediated by a spin flip, and E 1+ = S 1+ = . Dy- 
namically, quadrupole deformations may arise through 
the interaction of the photon with the pion cloud^^ 
or through the one-gluon exchange mechanism^ . Simul- 
taneous description of both R EM and R S m is achieved 
with dynamical models that include pion-nucleon inter- 
actions explicitly, supporting the idea that most of the 
quadrupole strength in the ATA (1232) transition is due 
to meson effects. At the real photon point over 30% of 
the transition strength is due to pion effects. 
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Fig. 1. Rem and Rsm extracted from exclusive re- 
actions p(e, ep)7T° using unitary isobar and fixed-t 
dispersion relations as analysis tools. 



At asymptotic momentum transfer, a model-independent 
prediction of helicity conservation requires Rem = 
E 1+ /M 1+ — > +1. An interpretation of Rem in terms of 
a quadrupole deformation can thus only be valid at low 
momentum transfer. Results of the multipole analysis 
of the JLab data 1 ^ 1 ^^^ as well as low Q 2 data from 
MAMl mm , Bates 111 and LEGS 1111 are shown in Figjl] 
A consistent picture emerges from these precise data: (1) 
The magnetic transition form factor drops much faster 
with Q 2 than the elastic form factor G v m which follows 
approximately the dipole form G D . (2) Rem remains 
negative, small and nearly constant in the entire range 



< Q 2 < 6 GeV 2 . (3) There are no indications that lead- 
ing pQCD contributions are important as they would re- 
sult in Rem —> +1 ? and (4) R S m remains negative, 
and its magnitude rises with Q 2 . 

Ultimately, we want to come to a QCD description 
of these important nucleon structure quantities. In re- 
cent years significant effort has been extended towards 
a Lattice QCD description of the N A transition^ 
Within the still large error bars, both quenched and un- 
quenched calculations at Q 2 < 1.5 GeV 2 with pion masses 
of 400 MeV are consistent with a small negative value of 
Rem in agreement with the data. For R SM there is a 
discrepancy at low Q 2 in both quenched and unquenched 
QCD calculation while the rise in magnitude of Rsm with 
Q 2 observed in the data is quantitatively reproduced in 
full QCD at Q 2 >lGeV 2 . 




Q 2 (GeV 2 ) 

Fig. 2. Magnetic transition form factor G M ex- 
tracted from exclusive reactions pion production, 
normalized to the dipole form. 



3 The Roper resonance Pn(1440) - a puz- 
zle resolved 

The 7V(1440)Pn resonance (aka "Roper") has been 
a focus of attention for the last decade, largely due to 
the inability of the standard constituent quark model to 
describe basic features such as the mass, photocouplings, 
and Q 2 evolution. This has led to alternate approaches 
where the state is treated as a gluonic excitation of the 
nucleon^, or has a small quark core with a large meson 
cloudim, or i s a hadronic molecule of a nucleon and a a 
meson^ . 

Given these different theoretical concept for the struc- 
ture of the state, the question "what is the nature of 
the Roper state?" has been a focus of the N* program 
with CLAS. The state couples to both Ntt and Ntttt final 
states. It is also a very wide resonance with about 350 
MeV total width. Therefore single and double pion elec- 
troproduction data covering a large range in the invariant 
mass W, with full center-of-mass angular coverage are key 
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in extracting the transition form factors in a large range 
of Q 2 . As an isospin I—\ state, the P 1 1 (1440) couples 
more strongly to n7r + than to p7r°. Also contributions of 
the high energy tail of the A (1232) are much reduced in 
that channel due to the I= \ nature of the A(1232). 




W(GeV) W(GeV) W(GeV) 

Fig. 3. W dependence of the three lowest Legendre 
moments from mr + angular distributions at fixed 
Q 2 = 2.05 GeV 2 . The dotted line indicates the 
cross section when the amplitudes of the Pi i (1440) 
are set equal 0. 

A large sample of differential cross sections and polarized 
beam asymmetry data from CLAS^ have been analyzed 
using the fixed-t dispersion relations approach (DR) and 
the unitary isobar model (UIM) 32 . Some of the features 
of the data may best be seen in the Legendre moments 
determined in fits to Response functions can be expressed 
in terms of Legendre polynomials, e.g. the azimuthal an- 
gle independent part of the differential cross section can 
be written as: 

oo 

(Jt + ea L = Dj +L Pt {cosGl ) . 

Figure [3] shows the lowest Legendre moments for this 
response functions. The transverse and longitudinal 
electro-coupling amplitudes A 1/2 and S 1/2 of the tran- 
sition to the N (1440) P 1 1 resonance are extracted from 
fits to the data 34 . They are shown in Fig. [JJ 




3 4 

Q 2 (GeV 2 ) 

Fig. 4. Transverse electrocoupling amplitude for 
the Roper Pn(1440) (left panel). The full cir- 
cles are the new CLAS results. The squares are 
previously published results of fits to CLAS data 
at low Q 2 . The right panel shows the longitudi- 
nal amplitude. The bold curves are all relativist ic 
light front quark model calculations The thin 
solid line is a non-relativistic quark model with a 



vector meson cloud 1 — , and the thin dashed line is 
for a gluonic excitation^. 

At the real photon point A 1/2 is negative. The CLAS 
results show a fast rise of the amplitude with Q 2 and a 
sign change near Q 2 = 0.5 GeV 2 . At Q 2 = 2GeV 2 the 
amplitude has about the same magnitude but opposite 
sign as at Q 2 = 0. It slowly falls off at high Q 2 . This 
remarkable behavior of a sign change with Q 2 has not 
been seen before for any nucleon. The longitudinal cou- 
pling S1/2 is smaller than the transverse one. The first 
results for the transition form factors of the Roper have 
recently been obtained in unquenched QCD 39 . The hy- 
brid baryon model is clearly ruled out for both ampli- 
tudes. At high Q 2 both amplitudes are qualitatively de- 
scribed by the light front quark models, which strongly 
suggests that the Roper is indeed a radial excitation of 
the nucleon. The low Q 2 behavior is not well described 
by the LF models and they fall short of describing the 
amplitude at the photon point. This indicates that im- 
portant contributions, e.g. meson-baryon interactions at 
large distances may be missing. In the light front system 
the Dirac and Pauli form factors Fl and F2 can be inter- 
preted as two-dimensional charge densities 36 . Figure^ 
shows the transition charge density in transverse impact 
parameter space. 
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Fig. 5. Transverse transition charge densities in the 
light front system for the proton to Roper transi- 
tion in the light-front helicity +1/2 -> +1/2. The 
light area indicates dominance of up quarks (pos- 
itive charge) , the dark area is due to down quarks 
(negative charge). 



4 The iV(1535)Sn state 

This state has been studied in the pr] channel where 
it appears as a rather isolated resonance near the Nrj 
threshold. Data from CLAS 1 ^ 1111 and Hall C 42 have 
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provided a consistent picture of the Q 2 evolution ob- 
tained from r] electroproduction data alone. There are 
two remaining significant uncertainties that need to be 
examined. The first one is due to the branching ratio of 
the 5n(1535) —>pr), the second one is due to the lack of 
precise information on the longitudinal coupling, which 
in the pr] channel is usually neglected. 

The PDG gives a range of (3^ G = 0.45 - 0.60. 
Since this state practically does not couple to channels 
other than Nrj and Nir, a measurement of the reaction 
ep — > en+n will reduce this uncertainty. Also, the Ntt 
final state is much more sensitive to the longitudinal 
amplitude due to a strong £ n — P n interference term 
present in the Ntt channel. Adjusting f3 N7T = 0.45 and 
Pnti — 0-45 brings the two data sets into excellent agree- 
ment as shown in Fig. [6j 



two amplitudes is given by: 




3 4 

Q 2 (GeV 2 ) 



Fig. 6. The transition amplitude A 1 / 2 for the 
5*11 (1535). The big full circles are from the 
analysis of the CLAS n7r + and data^^. 
The other data are from the analysis of pr] 
data Ell EH EU rp^ e theory curves are from var- 
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ious constituent quark models 



An important advantage of the Ntt channel is that it 
is also sensitive to the longitudinal transition amplitude 
which is due to a significant s—p wave interference with 
the nearby p-wave amplitude of the P n (1440). This can 
be seen in the multipole expansion of the Legendre mo- 
ment Dq T : 



D L n 



K 



Re(E 0+ SZ_+S 0+ M;_). 



5 Helicity structure of the L>i 3 (1520) 

A longstanding prediction ™ of the dynamical con- 
stituent quark model is the rapid helicity switch from 
the dominance of the A 3 / 2 at the real photon point to 
the dominance of the A 1 / 2 amplitude at Q 2 > 1 GeV 2 . 
In the simple non-relativistic harmonic oscillator model 
with spin and orbit flip amplitudes only, the ratio of the 



AD13 -| r y 2 

where a is a constant adjusted to reproduce the ratio at 
the photon point where A 1/2 is very small. 
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Q 2 (GeV 2 ) 

Fig. 7. Electro-coupling amplitudes A\/ 2 (left) and 
A 3/2 (middle) and S 1/2 (right) for the Di 3 (1520). 
Model curves as in Fig. 6. 

Figure [7] shows the CLAS results for the electro- 
couplings. The A 3/2 amplitude is large at the real photon 
point and decreasing rapidly in strength with increas- 
ing Q 2 . The A 1 / 2 amplitude increases rapidly in mag- 
nitude with increasing Q 2 , before falling off slowly at 
Q 2 > 1 GeV 2 . The rapid switch in the helicity structure 
can best be seen in the helicity asymmetry defined as 

Ahel = (^-1/2 " 



"^3/2)/ (^1 



/2 




A 2 /2 ), and shown in Fig. [8j 

+ I ( 



Q (GeV z ) 

Fig. 8. Helicity asymmetry for the photon induced 
proton-Di3 transition vs Q 2 . 



6 Search for new nucleon states 

The experimental program makes use of the CEBAF 
Large Acceptance Spectrometer (CLAS) ^ which pro- 
vides particle identification and momentum analysis in a 
polar angle range from 8° to 140°. The photon energy 
tagger provides energy-marked photons with an energy 
resolution of = 10~ 3 . Other equipment includes a 
coherent bremsstrahlung facility with a goniometer for 
diamond crystal positioning and angle control. The fa- 
cility has been used for linearly polarized photons with 
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polarizations up to 90%. There are two frozen spin po- 
larized targets, one based on butanol as target material 
(FROST), and one using HD as a target material (HD- 
Ice). The latter is currently under construction. FROST 
has already been operated successfully in longitudinal po- 
larization mode and will be used in transverse polariza- 
tion mode in 2010. The HD-Ice target will be used as 
polarized neutron (deuteron) target in scheduled runs in 
2010 and 2011. Circularly polarized photons can be gen- 
erated by scattering the highly polarized electron beam 
from an amorphous radiator. 

It is well understood that measurements of differential 
cross sections in photoproduction of single pseudoscalar 
mesons alone results in ambiguous solutions for the con- 
tributing resonant partial waves. The N* program at 
JLab is therefore aimed at complete, or nearly complete 
measurements for processes — > ttN, r]p, K + Y and 
jn — > 7r7V, KY . Complete information can be obtained 
by using a combination of linearly and circularly polar- 
ized photon beams, measurement of hyperon recoil polar- 
ization, and the use targets with longitudinal and trans- 
verse polarization. The reaction is fully described by 4 
complex parity conserving amplitudes, requiring 8 combi- 
nations of beam, target, and recoil polarization measure- 
ments for an unambiguous extraction of the production 
amplitudes. If all possible combinations are measured, 
16 observables can be extracted. In measurements that 
involve nucleons in the final state where the recoil polar- 
ization is not measured, 7 independent observables can 
be obtained directly, and the recoil polarization asym- 
metry P can be inferred from the double polarization 
asymmetry with linearly polarized beam and transverse 
target polarization. 
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Fig. 9. Integrated total cross section data for K + A 
and K + Y>° channels. The lines represent the 
fit results of the Bonn-Gatchina coupled-channel 
analysis. The dashed line shows the energy- 



dependence of the Pi 3 partial wave which indi- 
cates the presence of two Pi 3 resonances at 1720 
MeV and at 1900 MeV, respectively. In this anal- 
ysis, the new Pi3(1900) contributes a significant 
fraction to the total K + A cross section. 

6.1 Search for new states in strangeness chan- 
nels 

A large amount of cross section data have been 
collected in recent years on the KA and if E photo- 
production These data cover the nucleon res- 
onance region in fine steps of about 10 MeV in the 
hadronic mass W, and nearly the entire polar angle range. 
The integrated cross section shown in Fig. [9] reveals a 
strong bump around W = 1900 MeV which may be the 
result of s-channel resonance production. However, more 
definite conclusions can only be drawn when more po- 
larization observables are included in the analysis. Pre- 
liminary results are available from the beam asymmetry 
in A and E production with a linearly polarized pho- 
ton beam. A sample of the data is shown in Fig. 10 
Nearly full angle coverage with excellent statistics has 
been achieved. Beam asymmetry data are particularly 
important for identifying spin-parity assignment of con- 
tributing resonances. These data will span the resonance 
mass region from threshold up to W = 2.3 GeV. 

I Photon Asymmetry 1 .95GeV yp -> K*X° | 




Photon Asymmetry 1 95GeV 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos(6)Kcm 

Fig. 10. Sample of preliminary CLAS data for the 
beam asymmetry of the KY>° final state. 

In addition to precise KA and if E cross section data re- 
coil polarization and polarization transfer data have been 
measured The recoil polarization data in the K + A 
sector showed a highly unexpected behavior: The spin 
transfer from the circularly polarized photon to the A 
hyperon is complete, leaving the A hyperon 100% polar- 



ized, as can be seen in Fig. 11 The quadratic sum of all 
polarization components R 2 = P 2 +Cj;+C^, which has an 
upper bound of 1, is consistent with R 2 = 1 throughout 
the region covered by measurement. At first glance this 
result seems adverse to the idea that the KA final state 
has a significant component of N* resonance associated 
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with. However, the analysis of the combined CLAS dif- 
ferential cross section and polarization transfer data by 
the Bonn-Gatchina group 55 shows strong sensitivity to 
a P 13 (1900) candidate state. The decisive ingredient in 
this analysis are the CLAS spin transfer coefficients C x 
and C z . Note that the peak observed in the K + A data 
seen near 1900 MeV in Fig. [9] was originally attributed 
to a D 13 (1900) resonance before the spin transfer data 
were included in the fit. A P 13 (1900) is listed as a 2-star 
candidate state in the 2008 edition of the RPP m . If this 
assignment is confirmed in future analyses which should 
include additional polarization data, the existence of a 
7V(1900)Pi3 state would be strong evidence against the 
quark-diquark model but would support the symmet- 
ric SU(6) constituent quark model. 

First differential cross sections of the channel 7p — > 
covering the mass range W > 2100 MeV have 
been measured recently Several other strangeness 
channels such as -> K°A, K*°A, i^°E°, and 
K + £~(1385) are currently being studied to search for 
new states on neutron targets using circularly and lin- 
early polarized photon beams. For most of these chan- 
nels a complete kinematical reconstruction in the neutron 
rest frame is possible, thus eliminating the effect of Fermi 
motion in the deuteron nucleus. Recoil polarization mea- 
surements in all hyperon final states are also available for 
the analysis. 
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Fig. 11. Spin transfer from the polarized photon to 
the final state A. The spin transfer data were fit- 
ted by the Bonn-Gatchina group simultaneously 
with the CLAS K + A and K + Z differential cross 
section data. A Pi3(1900) state is required for the 
best fit to the cross section and the spin transfer 
data. 



7 Search for new cascade baryons 

The production of 2 hyperons, i.e. strangeness S = 
—2 excited states, presents another promising way of 
searching for new baryon states. The advantages of 
cascade hyperon states are due to the expected narrow 
widths of these states compared to S = 0, and S = — 1 
resonances. Possible production mechanisms include t- 
channel K or K* exchanges on proton targets with an 
excited hyperon Y* (A* or £*) as intermediate state 
and subsequent decays F* — > K + E* and 2* — > Ett or 
2 — > A(Y0)K. Missing mass technique may be used to 
search for new states in the reaction — » K + K + X if 
the state is sufficiently narrow to be observed as a peak 
in the missing mass spectrum. However, an analysis of 
the final state is needed to assign spin and parity to the 
state. Data from CLAS taken at 3.6 GeV beam energy 
show that one can identify the lowest two cascade states 
this way. 60 . To identify the higher mass states higher 
energy is needed. 

Another approach to isolate excited cascades and de- 
termine their spin-parity, is to measure H with additional 
particles in the final state. The invariant mass of the 



S°7r~ system is displayed in Fig. 12 and shows the first ex- 
ited state 2(1530) and indications of additional structure 
near 1620 MeV. A state near that mass is predicted ^ 
as a dynamically generated En system. The data have 
insufficient statistics and were taken at too low beam en- 
ergy to allow further investigations. New data taken at 
5.7 GeV electron energy and with higher statistics are 
currently being analyzed, and should allow more definite 
conclusions on a possible new 5 state at that mass. 



M:1.5392± 0.0014 
r:0.015± 0.005 
N: 71 ± 26 
X 2 /Ndf: 30.0/31.0 




1.5 1.55 1.6 1.65 1.7 1.75 1 



M(sV) (GeV/c 2 ) 



1.85 1.9 1.95 



Fig. 12. Invariant mass of the Ett~ final state. The 
2(1530) is clearly seen. A possible structure near 
1620 MeV could be related to a dynamically gen- 
erated 2 state 
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8 Strangeness electroproduction 

Significant effort has been devoted to the study of 
electroproduction of hyperons asa complementary 

means of searching for new excited nucleon states. Fig- 
ure [13] shows the dependence of the response functions 



on the hadronic mass W. 



P, E and double polarization asymmetries for beam- 
recoil polarization O y > , and target-recoil polarization T z > 
are shown. Other observables, e.g. target asymmetry T, 
beam-target asymmetries E, F, and G, H will be mea- 
sured as well. The projections are for one photon energy 
bin out of over 25 bins. The solid line is the projection 
of the kaonMAID code m . 
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Fig. 13. Separated response functions for K^Yf 
electroproduction. The response function a u and 
gtt show strong resonance- like structure near 
W — 1.9GeV, which is not reproduced by any 
model that include known resonances. 

The structure function a u shows a resonance-like en- 
hancement near 1.85 - 1.90 GeV in all angle bins. A 
comparison with model calculations reveals large discrep- 
ancies. None of the models that include known nucleon 
resonance couplings to i^E, is able to get the normal- 
ization correct. This leaves much room for yet to be 
identified resonance strength. 

9 Outlook 

The power of the polarization measurements on the 
extraction of resonant amplitudes will be brought to bear 
fully when double and triple polarization data are avail- 
able making use of longitudinally and transversely polar- 
ized proton and neutron targets. Data from CLAS on 
the polarized proton target (FROST) combined with lin- 
early and circularly polarized photons are currently in the 
analysis stage. Measurements with transverse polarized 
proton targets are planned for 2010. 

Measurements of polarization observables on polar- 
ized neutrons are planned for 2010/2011. Some projected 
data using the HD-Ice target facility in the CLAS detec- 



tor are shown in Fig. 14 Single polarization observables 



Fig. 14. Simulated data samples for the process 
771 — > K° A. The uncertainties are estimated with 
the expected run conditions for the upcoming po- 
larized target run using the HD-Ice facility at 
JLab. 



10 Conclusions 

With the recent precise data on pion and eta electro- 
production, combined with the large coverage in Q 2 , W, 
and center-of-mass angle, the study of nucleon resonance 
transitions has become an effective tool in the exploration 
of nucleon structure in the domain of strong QCD and 
confinement. We have learned that the A (1232) exhibits 
an oblate deformation. The multipole ratios R E m an 
R S m show no sign of approaching the predicted asymp- 
totic behavior, which provides a real challenge for model 
builders. Lattice QCD simulations do not yet reach the 
high Q 2 covered by the data. The latest data from 
CLAS on charged pion production reveal a sign change of 
the transverse amplitude for the N- Roper transition near 
Q 2 = 0.5 GeV 2 , and give strong evidence for this state as 
the first radial excitation of the nucleon. The hard transi- 
tion form factor of the fin (1535) previously observed only 
in the prj channel is confirmed in the nir + channel, which 
also allows us to extract the so far unmeasured longitu- 
dinal amplitude S 1 / 2 - The £) 13 (1520) clearly exhibits the 
helicity flip behavior predicted by the constituent quark 
model. In the light front frame the F 1 and F 2 transi- 
tion form factors that are obtained from the resonance 
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helicity amplitudes, can be interpreted as 2-dimensional 
light cone transition charge densities. This development 
marks real progress in our understanding of the spatial 
and spin structure of the nucleon in the domain of strong 
QCD. 

The search for new excited states of the nucleon has 
entered a new and exciting phase. The main goal of com- 
plete or even over determined measurements for pseudo- 
scalar meson production is now in reach with the avail- 
ability of energy-tagged linearly and circularly polarized 
photon beams and the construction of polarized proton 



and neutron targets that can be used in longitudinal and 
in transverse polarization modes. Moreover, for processes 
involving hyperons in the final state, the hyperon recoil 
polarization can be measured due to the large analyzing 
power of the hyperon weak decay A^pir~. 

Authored by The Southeastern Universities Research 
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AC05-84ER40150 . The U.S. Government retains a non- 
exclusive, paid-up, irrevocable, world-wide license to pub- 
lish or reproduce this manuscript for U.S. Government 
purposes. 
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